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Abstract: We have investigated the transmission and reflection properties 
i—i of structures incorporating left-handed materials with graded index of 

refraction. We present an exact analytical solution to Helmholtz' equation 

• i-H for a graded index profile changing according to a hyperbolic tangent 

O , function along the propagation direction. We derive expressions for the field 

intensity along the graded index structure, and we show excellent agreement 
c/3 between the analytical solution and the corresponding results obtained by 

accurate numerical simulations. Our model straightforwardly allows for 
c/3 arbitrary spectral dispersion. 
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1. Introduction 

During the last decade, a new class of artificial composite materials called electromagnetic 
metamaterials has emerged. Such metamaterials are structured at a subwavelength level and 
have electromagnetic resonators or "particles" such as split-ring resonators and nanowires as 
their structural units [ 1 1. Through the appropriate design of the electric and magnetic resonances 
of these subwavelength particles, it is possible to create materials with electromagnetic prop- 
erties that are unknown for conventional materials. Particularly, there has been considerable 
interest in metamaterials with simultaneously negative permittivity and negative permeability, 
which are called left-handed materials. It was shown in the theoretical work of Veselago [2| 
that left-handed materials exhibit a number of remarkable properties, including negative index 
of refraction (and, hence, negative phase velocity), inverse Doppler effect, and radiation ten- 
sion instead of pressure. All these properties stem from the fact that the Poynting vector in 
these materials is antiparallel to the wavevector, i.e., the electric field, the magnetic field and 
the wavevector of a plane electromagnetic wave form a left-handed system of reference. 

The concept of metamaterials was brought to practical implementation by a number of works 
by Pendry [ 3 , 4|, who suggested split-ring resonators and wire arrays as the first metamaterial 
particles, and the first experimental demonstration of a left-handed material was published in 
2001 [5|. Split-ring resonators and nanowires are still widely used in the microwave domain 
and are now well understood, but many other particles such as slab-wire pairs, fishnets 1 1 1 and 
coupled split-ring resonators [6| have been designed. One of the main challenges in the field of 
optical metamaterials is the reduction of the relatively high losses associated with their resonant 
behavior. The development of metamaterials with permittivity and/or permeability negative or 
smaller than unity has also led to a revolution in our understanding of optics; several remarkable 
applications have been proposed, including superlenses and hyperlenses that enable imaging 
below the diffraction limit [7, 8 9 1, waveguides that can stop light 1 10 1, miniaturized photonic 
devices such as Fabry-perot resonators and waveguides [11, 12 1, and even invisibility cloaks 
through the technique of transformation optics lfT3l[T4l . 

In this paper, we want to consider the transmission and reflection properties of graded index 
optical structures with a gradual transition from a right-handed to a left-handed material. Such 
structures were studied in the framework of metamaterial gradient index lenses by a few au- 
thors |[T5l[T6irT7l[T8ll . who have shown that this provides an additional degree of freedom that 
can be used, among others, to reduce geometrical aberrations; a gradient index metamaterial 



lens was also demonstrated experimentally by Smith 1191 . These works addressed the propaga- 
tion through graded index structures with geometrical optics and therefore neglect any reflec- 
tions that may occur on such structures. This issue has been taken up only very recently |20|. 

Here we present an exact analytical solution of Helmholtz' equation for the propagation 
of electromagnetic waves through a graded index metamaterial structure. We choose a graded 
index profile for which both the permittivity and the permeability vary according to a hyperbolic 
tangent function (see Fig. [TJ. In Sec. [2] we review the field equations that the electric and 
magnetic fields must satisfy in an inhomogeneous medium and we transform these equations 
to simple wave equations with an effective wavenumber. In Sec. [3] we present the analytical 
solution for the waves propagating in this structure. Finally, we compare the obtained analytical 
solution with results obtained from accurate numerical simulations based on a finite element 
method in Sec. [4] 
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Fig. 1. Propagation of an electromagnetic wave through a graded index structure with a 
hyperbolic tangent profile. This is the index profile assumed in this paper. 



2. Field equations 

We start our analysis from Maxwell's equations, and we search for fields that are periodic in 
time according to a exp(— icot ) dependency. Furthermore, we assume that the effective medium 
approximation can be made and that the materials are isotropic, so that their optical properties 
can be described by the effective dielectric permittivity and the effective magnetic permeability. 
For most metamaterials, the effective medium assumption is valid, because their constituents 
elements are on the subwavelength level. The geometry of the problem is illustrated in Fig.fTI 
The electric field is directed along the y-axis, E(r) = E(x)e y , whereas the magnetic field is 
directed along the z-axis, H(r) = H(x)e z . The propagation direction of the wave is along the 
x-axis. Since the fields depend only on the x-coordinate, we have 



dE 

dx 
dH 

dx 



= icofiH(x), 

= icoeE(x), 



(1) 
(2) 



where e = e(co,x) and ji = jj,(co,x) are the frequency-dependent dielectric permittivity and 
magnetic permeability, respectively. 

We can easily eliminate either the electric field or the magnetic field from Eqs. (fl]i-(J2|l ETI : 
this yields an ordinary differential equation for either E(x) or H(x): 



d 2 E 1 djl dE 
dx 2 jX dx dx 



(0 2 LLEE(x) = 0, 



(3) 



or 



d 2 H 1 de dH 
dx 2 £ dx dx 



(0 2 neH(x) = 0. 



(4) 



These equations describe the propagation of electromagnetic waves through a medium of which 
the constitutive parameters vary along the propagation direction. The spatial dependency of 
the functions e(x) and fi(x) may be completely arbitrary, even on space scales faster than the 
wavelength of the radiation, on the condition of course that the effective medium approximation 
remains valid. 

The standard approach to the solution of Eqs. ([3])-(|4]l is to eliminate the first order terms 
by introducing the functions F(x) and G(x) instead of the functions E(x) and H(x) using the 
following transformations: 



E(x) = y/jip}F(x), H(x) 



In this way, we obtain the following wave equations for the functions F(x) and G(x); 
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These two equations can also be written as wave equations 
d 2 F .,,,„,, „ d 2 G 



G(x) = 0. 
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are the space-dependent effective wave vectors for electric and magnetic fields respectively. In 
case of the hyperbolic tangent profile for the functions e(x) and fi(x), Eqs. d6]l-(|7]l are generally 
reduced to the hypergeometric equation, allowing for analytical solution in terms of suitable 
hypergeometric functions. 

3. Analytical solutions of the field equations 

In this paper, we consider an inhomogeneous medium for which the effective permittivity and 
permeability vary according to a hyperbolic tangent function. We choose the hyperbolic tangent 
function here because we will show that it is possible to construct analytical solutions for this 
case. Furthermore, it is a convenient function that naturally provides asymptotic values of the 
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Fig. 2. We assume in this paper that the effective permittivity E e ff and permeability jUg// 
vary along the propagation direction according to a hyperbolic tangent function. 



constitutive parameters in the backward and forward directions and allows a detailed study of 
the limit of the abrupt transition as well. We use the antisymmetric functions (see Fig. [2]i 



\i = -non e ff((o)tanh(pi 



-eoe £ .//(a))tanh(px), 



(11) 



where p is a positive real parameter describing the steepness of the transition from the right- 
handed material at the left-hand side of the plane x — to the left-handed material at the 
right-hand side of the plane x = 0. There is no restriction on the functions ji e ff((o) and 
E e ff{(0) (except of course for such restrictions as the Kramers-Kronig relationships), so that 
our method allows for arbitrary spectral dispersion. The reader should note that the impedance 
Z = ZqZ(co) = J HofJ.eff((o) I £o£eff((o) is constant throughout the entire structure; as a result, 
there is no reflection on the graded interface between the two materials. 

We found that, for this particular graded index structure, the two differential equations Q-Q 
have a remarkably simple set of two independent exact solutions given by 



E(x)=E [cosh(px)] ±i P , 
where Eq and Hq are constant amplitudes, and 



H(x) = Hq [cosh (px] 



i±/l 



C0 2 £ eff (G))ll e ff((O). 



(12) 



(13) 



We repeat that these exact solutions are valid for arbitrary steepness p . Let us now choose the 

(14) 



solution with the minus sign in the exponent of the expression ( 12 1, i.e., 

E(x)=E [cosh(px)]~ i P . 



From the asymptotic expansion of the field E(x) in the limits x — > =F°°, 

E(T°°)=E e i p la2 e± iKX , 



(15) 



we see that the wave described by Eq. ([14]) is an electromagnetic wave with wavevector 
k_oo = + K"e x in the right-handed material far from the interface (x —> — °°). This is a wave 
that propagates in the +x direction, i.e., a wave propagating to the right. On the other hand, 
for x — s- +°°, the wave has wavevector k +0 o = — Ke x ; this represents a wave of which the phase 
fronts propagate in the — x direction. However, since we have a left-handed material for x > 0, 
the energy flux (Poynting's vector) is still propagating from left to right. This is perfectly con- 
sistent with the fact that there is no reflection on this structure. This is also apparent from the 
fact that | cosh(px)~' K > p \ = 1, so that \E(x)\ is constant throughout the structure. 



(a) analytical solution: 
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Fig. 3. Comparison of the analytical and numerical results for the electric field. We plot 
E(x) at t = 0. (a) Analytical solution for p = 10 um. (b) Numerical solution for p = 10 |i.m. 
(c) Analytical solution for p = 1 ixm. (b) Numerical solution for p = 1 |j.m. 



4. Comparison with numerical results 

In order to validate our exact analytical solution [Eq. (12)], we compare the waveforms with re- 
sults obtained from a direct simulation of Maxwell's equations. We use a finite element method 
(COMSOL Multiphysics) to discretize Maxwell's equations, using perfectly matched layers 
at the left and the right of the structure to close the simulation domain. For the simulation 
results shown in Figs. I5lb) and (d), we have used the following parameters: Aq = 1 (o.m and 
£ e ff (Ao) = jU e ff(Ao) = 1. Figures [3ta)-(b) and Figs.[3tc)-(d) are for different transition steepness. 
We see that there is excellent agreement between the analytical and numerical results. 

5. Conclusion 

We have investigated electromagnetic wave propagation through a graded index interface be- 
tween a right-handed and a left-handed material. We derived an exact analytic wave solution 
when the index of refraction varies according to a hyperbolic tangent function. This solution is 
valid for arbitrary steepness of the index transition, even when the effective constitutive param- 
eters vary on the scale of the vacuum wavelength where the traditional approximate methods 
(e.g., SVEA) cease to work. We have validated our analytical solutions by accurate numerical 
simulations using a finite element method. Our analytical model allows for arbitrary dispersion. 
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